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This  summary  of  Late-Time  Containment  Research  describes  the  various  tasks 
undertaken  during  1981.  Cavity  conditions  were  estimated  for  Red  Hot  as  a 
function  of  time,  and  consideration  of  connecting  tunnel  conditions  suggested 
that'  venting  from  the  cavity  was  energetically  rather  minor.  Olablo  Hawk 
diagnostic  data  permitted  the  Identification  of  the  second  LOS  flow  pulse  as 
steam  at  a  temperature  near  3500° K.  The  LOS  flow  In  Miners  Iron  was  found  to 
have  a  pressure  three  times  higher  than  on  Diablo  Hawk.  LOS  flow  studies  (over) 
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suggested  that  aerodynamic  forces  turn  the  helical  insert  into  the  flow  and 
thereby  Intercept  penetrating  flow  components.  Particle  velocity  records 
from  grout  spheres  tests  were  analyzed  to  show  that  the  conventional 
material  model  for  2C4  grout  is  unacceptable.  Probably  strain  rate 
effects,  material  damage  and  a  relaxation  time  will  be  required  to  give 
agreement  between  calculations  and  the  measurements.  An  algorithm  was 
developed  for  use  in  the  LANL  KRAK  code  to  improve  calculations  of  crack 
propagation  from  a  cylindrical  borehole.  This  specified  how  pressure  in 
the  borehole  should  be  treated.  Several  other  tasks  are  discussed  briefly. 
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SECTION  1 
INTRODUCTION 


This  Is  the  final  report  of  activities  undertaken  during 
calendar  year  1981  In  the  continuing  support  provided  by  S-CUBED  to 
the  Defense  Nuclear  Agency  In  the  general  area  of  Late  Time  Containment 
research.  This  work  was  supported  by  Contract  No.  0NAG01-81-C-0072 . 

As  In  previous  years  an  effort  has  been  made  during  1981  to 
report  technical  activities  In  a  timely  way  In  topical  reports  or  In 
letters  to  Interested  parties.  The  presentation  found  here  will,  where 
appropriate,  rely  heavily  upon  these  more  complete  discussions. 

Contributions  have  been  made  In  a  number  of  areas.  These 
include  the  analysis  of  past  events  and  suggestions  relevant  to  the 
design  of  future  experiments,  studies  of  line  of  sight  of  flow  condi¬ 
tions,  support  for  the  LS  experimental  series  conducted  by  Physics 
International,  support  for  the  grout  spheres  program  conducted  by  SRI 
International  and  a  number  of  miscellaneous  activities.  These  Items 
will  be  discussed  In  the  following  sections. 
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SECTION  2 

TEST  RELATED  ACTIVITIES 


Obviously  a  major  responsibility  of  the  Late  Time  Containment 
program  Is  to  provide  design  Information,  comment,  and  evaluation 
concerning  specific  nuclear  tests  which  have  been  conducted  or  which 
are  planned  by  DNA.  Contributions  In  several  such  areas  will  be 
discussed  below. 

2.1  RED  HOT/MISTY  JADE 

A  major  current  effort  In  the  DNA  containment  program  Involves 
the  development  of  a  containment  scenario  for  a  proposed  cavity  shot. 
Misty  Jade,  which  Is  In  many  respects  similar  to  Red  Hot.  As  a  result, 
an  extensive  program  Is  underway  to  understand  the  Red  Hot  experience. 
It  Is  generally  believed  that  this  understanding  Is  a  prerequisite  to 
the  development  of  a  satisfactory  containment  prospectus  for  Misty 
Jade.  The  status  of  understanding  of  Red  Hot  phenomenology  as  of 
May,  1981,  Is  summarized  In  a  classified  report  entitled  "Pressure 
Evolution  Following  a  Nuclear  Explosion  In  a  Cavity".  The  Introduc¬ 
tion,  Summary  and  Conclusion  section  of  that  report  are  reproduced 
below. 


"A  nuclear  explosion  In  a  pre-mi ned  hemispherical  cavity, 
such  as  the  Red  Hot  event  and  the  proposed,  very  similar  Misty  Jade 
event.  Is  of  Interest  for  the  measurement  of  cratering  effects. 
However,  such  an  explosion  Is  largely  decoupled  from  the  surrounding 
medium.  As  a  result  a  large  fraction  of  the  yield  remains  In  the 
cavity.  This  raises  concerns  about  the  satisfactory  containment  of 
the  high  pressure  cavity  gases. 

"This  report  analyzes  the  evolution  of  cavity  pressure  as 
more  and  more  of  the  surrounding  rock  is  equilibrated  with  the  cavity 
gas.  This  material  may  be  mixed  with  the  cavity  gases,  for  example, 
as  a  result  of  ablation  of  crater  ejecta  during  Its  flight  through 
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the  cavity  or  spallation  of  the  cavity  wall  as  a  result  of  heat  transfer 
from  the  cavity  gases.  The  other  major  source  of  material  for  cavity 
quenching  Is  chimney  formation. 

"The  cavity  pressures  at  various  states  of  mass  addition  for  Red 
Hot  are  estimated  in  Section  2.  It  is  concluded  that  tens  of  thousands 
of  tons  of  rock  are  required  to  get  the  cavity  pressure  down  to  safe 
levels  (below  overburden)  for  the  yield  and  cavity  volume  considered. 
There  Is  an  intermediate  regime  where  the  pressure  Is  highly  "buffered" 
changing  only  moderately  as  the  added  rock  mass  Increases  from  a  few 
hundred  tons  to  about  ten  thousand  tons.  The  buffered  level  Is  estima¬ 
ted  to  be  about  140  to  200  bars  for  Red  Hot,  assuming  20%  moisture.  The 
level  Is  estimated  at  80  to  110  bars  for  10%  moisture,  suggesting  a 
significant  advantage  In  selecting  a  dry  formation  for  any  future  cavity 
shot. 


"The  rates  of  removal  of  material  from  ejecta  and  from  the  cavity 
wall  are  considered  In  Section  3.  It  Is  concluded  that  ejecta  ablation 
could  not  reduce  the  Red  Hot  cavity  pressure  to  safe  levels  (below 
overburden)  and  that  cavity  wall  spallation  could  do  so  only  after  a 
number  of  hours.  The  pressure  is  expected  to  be  reduced  to  the  buffered 
level  within  several  seconds  to  several  minutes. 

"An  Inventory  was  made  of  the  possible  mechanisms  whereby  the 
cavity  pressure  may  have  been  relieved  In  Red  Hot.  Although  only 
limited  Information  has  been  obtained  by  reentry  work  to  date.  It  has 
been  tentatively  concluded  that  the  energy  vented  Into  the  access  tunnel 
and  Deep  Well  was  Insufficient  to  reduce  cavity  pressure  significantly. 
The  pressure  may  have  been  relieved  by  a  crack  opening  Into  a  porous 
medium  above  the  cavity  or  by  prompt  chimney  formation,  or  it  may  have 
been  sufficiently  reduced,  by  a  combination  of  prompt  ablation  and  the 
effect  of  relatively  dry  rock  adjacent  to  the  cavity,  to  enable  the 
residual  stress  field  to  contain  It.  It  Is  even  possible  that  the 
residual  stress  field  was  sufficient  to  contain  the  pressure  of  perhaps 
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400  bars  or  so  to  which  the  cavity  can  rapidly  drop  by  equilibrating  the 
exposed  sand  and  st?e1  known  to  be  present.  This  analysis  Is  presented 
In  Section  4. 

"The  Implications  of  the  Red  Hot  observations  and  the  analysis 
contained  In  this  report  for  Misty  Jade  are  also  discussed  In  Section  4. 
It  Is  concluded  that  containment  success  for  Misty  Jade  cannot  be  confi¬ 
dently  predicted  at  this  time.  Containment  prospects  can  be  maximized  by 
emplacing  It  In  the  driest  medium  possible  without  compromising  Its 
ability  to  produce  a  residual  stress  field. 

"It  should  be  stressed  that  this  Investigation  has  been  carried  out 
prior  to  a  planned  Red  Hot  reentry  and  Is  Intended  In  part  to  provide 
guidance  for  that  work.  The  findings  of  the  reentry  should  enable  us  to 
narrow  down  the  Red  Hot  scenario  much  more  than  Is  currently  possible  and 
to  make  predictions  for  Misty  Jade  with  greater  confidence." 

In  addition,  a  number  of  meetings  and  conversations  have  occurred 
In  which  a  number  of  recommendations  for  the  development  of  a  productive 
reentry  program  have  been  presented  to  Field  Command  representatives.  We 
continue  to  stand  ready  to  assist  with  the  analysis  and  Interpretation  of 
reentry  observations  as  they  become  available. 

2.2  DIABLO  HAWK  CONTAINMENT  SUWARY  REPORT 

The  experimental  data  related  to  LOS  flow  conditions  and  late  time 
gas  pressure  and  composition  have  been  collected,  analyzed  and  Inter¬ 
preted  In  a  classified  Containment  Summary  Report  for  the  Diablo  Hawk 
event.  Specifically  this  report  describes  the  LOS  plasma  flow,  the 
stemming  motion,  tunnel  environment,  residual  stress  and  stress  relaxa¬ 
tion  and  the  performance  of  associated  tunnel  systems. 

The  reverse  cone  used  for  Diablo  Hawk  was  thinner  than  that  used 
on  Mighty  Epic,  and  the  Iron  vacuum  vessel  was  outside  of  the  tuballoy 
segments  rather  than  Inside  of  them.  This  design  change  was  based  on 


calculatlonal  studies  and  problems  which  arose  In  the  construction  of  the 
Nighty  Epic  system. 

Many  diagnostic  measurements  were  made  to  determine  the  plasma 
characteristics  In  the  LOS  for  this  event.  In  general,  good  results  were 
obtained,  and  the  data  provide  the  basis  for  an  Improved  understanding  of 
the  flow  phenomenology.  Specifically,  sllfer,  fluid  coupled  plate,  and 
ported  pressure  systems  were  Installed,  and  good  data  were  obtained  from 
most  of  them. 

Relatively  weak  HMTW  cable  was  used  for  slifer  systems  on  Diablo 
Hawk.  This  cable  has  an  estimated  crush  strength  approximately  one-half 
that  of  the  standard  cables  used  on  Mighty  Epic  and  previous  events. 

These  sllfer  systems  clearly  demonstrate  that  the  plasma  conditions  on 
Diablo  Hawk  were  more  benign  than  In  any  previous  experiment.  The 
Brownlee  knee  was  seen  at  a  range  of  28  m,  and  the  data  show  strong 
evidence  of  a  second  plasma  pulse  arriving  somewhat  later.  Qualita¬ 
tively  speaking, the  second  pulse  was  more  promlnant  In  Diablo  Hawk  than 
It  was  In  Mighty  Epic,  but  It  Is  clearly  evident  In  the  sllfer  records 
from  both  events. 

Fluid  coupled  plate  and  ported  pressure  gauges  were  Installed 
between  '0  and  110  m  along  the  LOS  pipe.  All  gauge  systems  provided  data 
and  confirm  the  two  wave  character  of  the  plasma  flow  Indicted  by  the 
sllfers.  The  fast  pulse  has  an  apparent  velocity  of  1.3  cm/us;  the 
slower  pulse  a  value  of  0.42  cm/iis.  In  both  cases  the  velocities  are 
constant  with  range. 

This  diagnostic  Information  was  used  to  develop  a  flow  scenario  for 
the  second  pulse.  It  appears  to  have  an  Incident  pressure  of  approximate¬ 
ly  0.2  MPA,  a  density  of  1.2  x  10"*  gm,  a  sound  speed  of  0.15  cm/us, 
and  a  temperature  of  3500*  K.  The  gas  Is  thought  to  be  largely  steam 
generated  from  the  collapsing  grout  near  the  end  of  the  tubal loy  reverse 


cone. 


The  calculated  cavity  formation  and  grout  extrusion  process  sug¬ 
gested  that  stemming  materials  should  have  reached  approximately  95  m 
from  the  working  point.  On  the  other  hand,  on  reentry  it  was  found  that 
the  OAC  doors  were  completely  shattered.  Fragments  were  found  throughout 
the  region  slightly  before  the  DAC  up  to  the  MAC.  A  significant  quantity 
of  grout  and  three  short  sections  of  LOS  pipe  passed  through  the  OAC 
opening.  The  DAC  door  fragments  were  heavily  eroded  by  what  appears  to 
have  been  an  extreme  thermal  environment.  Fragments  appeared  to  have 
been  well  "fried". 

The  relative  motion  of  the  grout  at  the  LOS  and  at  the  tunnel 
walls  was  dramatically  Illustrated  by  the  observation  of  about  20  m  of 
displacement  of  LOS  pipe  sections  near  the  DAC  while  essentially  no 
relative  motion  occurred  between  tuff  and  grout  at  the  tunnel  wall  In 
the  same  region. 

Finally,  chemical  analysis  of  samples  from  the  DAC/MAC  area  showed 
that  fission  fragments  from  uranium  as  well  as  many  construction  mater¬ 
ials  were  present.  The  uranium  came  from  the  predicted  jetting  of  the 
tubal loy  cone,  the  fission  debris  was  expected  because  of  the  neutron 
flux  on  the  cone.  The  relative  amounts  of  various  isotopes  observed 
were  not  characteristic  of  bomb  debris.  There  was  no  evidence  of  any 
leakage  from  the  cavity  to  the  DAC  region. 

Pressure  gauges  had  been  Installed  In  the  doorwells  of  the  DAC  and 
MAC.  Both  provided  Information  which  indicates  that  the  DAC  doors 
failed  at  approximately  400  ms.  It  Is  speculated  that  the  stagnated 
gases,  primarily  from  the  second  pulse,  were  compressed  by  the  oncoming 
grout  column,  and  these  gases  provided  the  heat  transfer  medium  which 
caused  the  observed  erosion  upon  subsequent  door  failure. 

Perhaps  the  observed  grout  extrusion  15  m  beyond  the  point  at 
which  calculations  had  Indicated  extrusion  would  stop  was  related  to  the 
relatively  large  tunnel  used  In  the  Diablo  Hawk  compared  to  that  of 


Nighty  Epic.  True,  the  Invert  had  been  filled  with  tuff  matching  grout 
in  an  effort  to  make  the  configurations  appear  Identical.  Perhaps  tuff 
matching  grout  does  not  behave  quite  like  tuff  under  dynamic  loads. 

The  late  time  pressure  In  the  TAPS  region  reached  a  pressure  of  24 
psla  2.5  hours  after  Diablo  Hawk,  and  it  had  decayed  to  ambient  at  the 
end  of  six  days.  In  addition,  gas  chromatographic  measurements  were 
made  on  samples  collected  from  behind  the  TAPS  which  indicated  that  the 
hydrogen  concentration  was  65%  ,  CO  concentration  was  1%  ,  oxygen  and 
nitrogen  totalled  approximately  30%.  Chemical  arguments  were  presented 
based  on  the  reactions  of  Iron  and  water  and  iron  and  C02  which 
suggest  that  the  temperature  behind  the  DAC  reached  4000-5000  K  as  gas 
was  compressed  by  the  oncoming  grout  column.  Reactions  with  iron  would 
have  changed  half  of  the  water  into  hydrogen  and  simultaneously  reacted 
most  of  the  C02  present  to  CO. 

Several  conclusions  can  be  drawn  from  the  study  of  the  late  time 
tunnel  environment.  The  first  Is  that  the  observed  explosive  mixtures 
are  normal  for  such  events  and  they  come  quite  naturally  from  the  shock 
closure  process,  LOS  pipe  wall  ablation,  and  known  chemical  reactions 
under  reasonable  thermodynamic  conditions.  The  most  Important  point  is 
that  these  explosive  gas  components  arise  In  processes  occurring  well 
removed  from  the  cavity.  The  observations  of  explosive  gases  beyond  the 
stemming  of  HLOS  nuclear  events  at  late  time  does  not  Imply  any  cavity 
leakage  Into  the  tunnel  system. 

An  effort  was  made  to  measure  the  residual  stress  through  a 
hydrofracture  system  installed  between  53  and  84  m  of  the  working 
point.  Only  the  system  at  68.5  m  was  not  heavily  damaged  during  the 
event.  It  clearly  shows  a  significant  increase  In  fracture  pressure 
within  5  minutes  of  the  explosion.  Subsequently  It  too  failed.  Nine 
days  after  the  event  the  system  damage  was  repaired,  and  hydrofracture 
measurements  were  then  taken  on  all  three  systems.  No  significant 
residual  stress  was  observed  at  that  time.  It  is  apparently  true  that  a 
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significant  residual  stress  was  developed  at  the  time  of  the  explosion, 
but  that  some  relaxation  process  had  eliminated  this  stress  within  a 
nine  day  period. 

2.3  MINERS  IRON  DIAGNOSTIC  REVIEW 

The  Miners  Iron  diagnostic  data  available  In  the  spring  of  1981 
was  analyzed  and  compared  with  Diablo  Hawk.  This  comparison  Is  part  of 
the  effort  required  for  the  preparation  of  the  Miners  Iron  containment 
summary  report  which  will  be  written  when  the  reentry  has  been  completed. 
Several  factors  were  evident  In  this  comparison.  The  sllfers  again 
showed  a  clear  Brownlee  knee  but  It  was  at  a  range  of  40  m  compared  to 
28  m  for  Diablo  Hawk.  It  was  also  apparent  that  muffler  Interactions 
with  the  first  plasma  pulse  caused  sllfer  signals  which  had  not  been 
seen  heretofore.  Fluid  coupled  plate  and  ported  pressure  gauges  again 
provided  particularly  good  data.  Two  pulses  travelling  with  velocities 
quite  similar  to  those  of  Diablo  Hawk  were  again  observed.  The  pressure 
amplitudes,  however,  were  approximately  a  factor  of  three  larger  than 
those  seen  on  Diablo  Hawk,  but  they  were  still  slightly  below  the 
average  experience  as  determined  by  Sandla  measurements. 

In  this  case  the  expectation  was  that  grout  extrusion  would 
destroy  the  closest  closure  at  98  m.  Reentry  shows  that  door  survived. 

Unexpectedly  high  pressures  were  observed  In  front  of  and  between 
the  closures  and  In  front  of  the  TAPS  within  a  few  minutes  after  Miners 
Iron.  Pressures  as  high  as  300-500  psl  are  estimated  to  have  existed 
throughout  these  regions  within  the  first  hour  after  the  event.  Unfortu¬ 
nately  no  late  time  pressure  or  temperature  Information  Is  available, 
and  no  gas  quality  measurements  comparable  to  those  from  Diablo  Hawk  and 
previous  events  are  available  from  Miners  Iron.  As  a  result,  It  will  be 
Impossible  to  discuss  the  chemical  Implications  of  the  late  time  pressure 
environment  In  the  MAC/TAPS  region  on  Miners  Iron  In  satisfactory  detail. 


2.4  HURON  UNDING  EXPERIMENTAL  RECOWENDATIONS 


A  number  of  recommendations  were  made  for  diagnostic  systems  to  be 
Included  In  the  Huron  Landing  event.  Specifically,  It  was  suggested  that 
sllfers  employ  HMTW  cable  as  In  Miners  Iron  and  Diablo  Hawk  and  that  at 
least  one  sllfer  be  Installed  quite  close  to  the  LOS  pipe  to  better 
Indicate  flow  conditions  within  the  LOS. 

Huron  Landing  Is  the  first  HLOS  event  to  Include  a  helical  insert 
to  attenuate  plasma  and  penetrating  debris.  It  was  recommended  that 
particular  efforts  be  made  to  Install  fluid  coupled  plate  systems  near 
45  m  to  diagnose  the  performance  of  this  new  feature.  Gauges  should  be 
located  between  and  under  the  segments  themselves.  Special  care  should 
be  taken  to  avoid  electrical  noise  and  to  allow  for  relatively  non- 
uniform  pipe  expansion. 

It  would  be  well  to  devise  some  kind  of  door  Impact  monitor.  If 
possible,  and  to  Install  MAC  doorwell  pressure  monitors  of  the  sort  used 
on  Diablo  Hawk.  These  might  provide  interesting  data  about  the  late  time 
LOS  pressure  and  the  time  of  failure  of  the  doors  if  grout  Impact  were  to 
occur. 


Finally,  special  efforts  were  recommended  to  obtain  good  gas 
quality  measurements  and  late  time  pressure  and  temperature  indications 
for  Huron  Landing.  The  Diablo  Hawk  data  suggested  the  chemical  mechan¬ 
ism  In  which  Hj/HgO  and  CO/CO 2  ratios  are  particularly  significant. 

He  hope  that  gas  samples  can  be  taken  and  subsequently  analyzed  to  pro¬ 
vide  comparable  data  from  Huron  Landing. 

2.5  ADDITIONAL  SUGGESTIONS 

Additional  rather  short  and  Informal  analyses  were  made  for  the 
Huron  Landing/Diamond  Ace  configuration  which  were  reported  to  the  COR 
through  correspondence.  Specifically,  investigations  were  made  of  the 
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posslbllty  of  using  high  explosive  detonation  products  In  the  vicinity  of 
the  muffler  to  slow  down  the  first  plasma  pulse  enough  to  Insure  Its 
capture  by  the  DAC  systems.  It  appeared  that  this  could  be  accomplished 
If  approximately  19  kg  of  explosive  were  detonated  In  the  muffler  region. 

Consideration  was  given  to  the  possibility  of  constructing  the 
Diamond  Ace  LOS  pipe  from  nested  segments  of  cylindrical  pipe,  each 
segment  of  which  would  be  a  stock  pipe  size.  An  apparently  attractive 
scheme  for  constructing  a  pipe  and  providing  helical  inserts  was  proposed. 


SECTION  3 

ANALYSIS  OF  JETTING  SIMULATION  EXPERIMENTS 


This  section  presents  a  brief  abstract  of  work  which  has  already 
appeared  In  a  detailed  topical  report  (Ref.  2)  prepared  during  the 
contract  period  as  well  as  some  additional  results  and  Insights 
obtained  after  the  release  of  that  report. 

3.1  REVIEW  OF  WORK  APPEARING  IN  TOPICAL  REPORT 

The  report  (Ref.  1 )  describes  calculations!  and  analytical 
studies  of  high  explosive  jetting  experiments  designed  to  simulate 
energy  Injection  into  a  nuclear  llne-of-slght.  Numerous  experiments 
were  performed  by  Physics  International  Co.  (Ref.  2)  using  Identical 
spherical  HE  sources  embedded  In  wet  sand.  Damage  was  assessed  by 
comparing  the  volume  and  depth  of  holes  produced  In  aluminum  target 
plates  at  the  end  of  each  pipe.  An  evacuated  steel  pipe,  19  mm  o.d. 
and  0.3  mm  thick,  was  the  standard  configuration.  Many  variations 
were  used. 

Axially  symmetric  experiments  produced  some  unexpected  results, 
for  example:  use  of  a  2  n  jacket  of  lead  around  a  standard  pipe 
produced  holes  about  twice  as  great  In  depth  and  volume  as  the  average 
of  the  standard  configurations;  use  of  an  air-filled  pipe  with  wall 
thickness  of  0.7  mm  produced  triple  the  standard  hole  size.  Two 
dimensional  calculations  for  the  Standard  and  Lead  Wrap  cases  which 
did  not  Include  the  thin  pipe  wall  were  totally  Inconsistent  with 
these  results.  Calculations  Including  the  wall  were  quite  consistent 
with  the  data  If  relative  calculated  energy  In  the  particulate  jet  of 
wall  material  formed  on  collapse  Is  compared  with  relative  measured 
hole  volume.  In  terms  of  absolute  comparison,  the  calculated  results 
are  consistent  only  If  the  cratering  efficiency  of  the  shower  of 
particulate  matter  Is  assumed  to  be  three  to  five  times  smaller  than 
that  of  a  single  particle  at  the  same  velocity  -  which  Is  considered 
reasonable. 
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Simplified  calculations  including  the  wall  were  performed  for  all 
three  axisymmetric  configurations  in  the  frame  of  reference  of  a  steady 
shock  to  provide  guidance  for  the  calculations  of  the  actual  configurations 
with  a  decaying  shock.  These  produced  the  unexpected  result  that  all 
three  jet  particulate  mass  and  energy  at  equal  rates  as  a  steady  state  is 
approached.  This  result  is  explained  by  an  adaptation  of  impulsive  jetting 
theory.  The  steady  state  results  put  constraints  on  explanations  of  the 
laboratory  results.  The  stronger  jetting  in  the  Lead  Wrap  case,  compared 
with  the  Standard,  is  attributed  to  slower  response  of  the  region  of  the 
source  of  the  jet  to  decay  in  the  driving  pressure  in  the  wet  sand  -  the 
Lead  Wrap  configuration  is  better  tamped.  The  0.7  mm  air  filled  pipe 
result  is  explained  analogously. 

Configurations  employing  helical  spiral  asymmetries,  0.3  to  0.5  mm 
thick,  also  produced  unexpected  results.  Spirals  resting  against  the 
outer  pipe  wall  resulted  in  essentially  no  change  in  target  damage,  whereas 
spirals  resting  against  the  inner  wall  virtually  eliminated  target  damage. 
As  a  consequence,  original  intuition  that  the  spiral  would  prevent  crater¬ 
ing  because  of  inertial  effects  was  abandoned.  A  postulated  mechanism 
attributing  the  spiral  effect  to  turning  of  the  unattached  spiral  into  the 
pipe  by  the  action  of  the  faster-arriving  vapor  flow  (generated  by  ener¬ 
getic  pipe  implosion  near  the  source)  is  presented  in  detail.  This 
mechanism  is  shown  to  be  consistent  with  all  data  acquired  so  far  in  the 
series  of  experiments.  (Note  added  in  proof:  So-called  scoping  tests 
fired  after  this  report  was  written  do  not  show  the  postulated  turning  of 
the  helix  into  the  flow.  Rather,  flash  X-ray  photographs  seem  to  show  a 
catastrophic  interaction  between  the  helix  and  unobservable  particulates 
which  is  accompanied  by  significant  conversion  of  kinetic  into  internal 
energy. ) 

3.2  FURTHER  ANALYSIS  OF  JETTING  PHENOMENA 

The  scenario  emerging  from  the  work  in  Ref.  1  envisions 
shock-driven  pipe  collapse  as  the  source  of  a  jet  of  condensate  or 
particulate  material.  This  jet,  which  comes  from  the  steel  liners, 
forms  the  observed  hole  in  the  target  plate  by  a  process  akin  to 
hypervelocity  Impact  cratering.  The  jet  consists  of  a  large  number  of 
small  particles  or  droplets;  interference  between  these  causes  the 
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"cratering  efficiency"  to  be  several  times  smaller  than  that  observed 
with  single  incident  hypervelocity  projectiles. 

Some  of  the  collapsing  liner  Is  vaporized,  especially  at  the 
higher  driving  pressures  near  the  source.  This  vapor  is  thought  to 
contribute  little  to  the  target  damage.  As  noted  above.  It  Is  also 
thought  (by  us  at  least)  to  be  responsible  for  the  turning  of  Internal 
spiral  helices  Into  the  path  of  the  particulate  jet,  eliminating 
target  damage  almost  entirely. 

More  recent  efforts  have  given  Increased  support  to  the  above 
scenario.  These  include  a  quasl-analytlc  study  of  phase  change  of  the 
collapsing  liners  and  a  review  of  the  Implications  of  data  from  the 
more  recent  LS-5  experiment  by  Physics  International  (PI),  Ref.  3d. 

We  now  describe  these. 

The  collapse  of  a  liner  (pipe  wall)  driven  by  a  shock  propaga¬ 
ting  parallel  to  the  pipe  axis  In  the  surrounding  medium  has  been 
treated  by  an  approximate  theoretical  model  presented  In  Appendix  A  of 
Ref.  1.  The  model  Is  an  adaptation  of  the  theory  of  Blrkhoff,  et  al 
(Refs.  3,4)  which  describes  steady  state  jetting  of  a  liner  driven 
Impulsively  by  a  layer  of  high  explosive.  The  present  analysis 
considers  a  liner  collapse  driven  by  the  pressure  field  behind  a 
constant-velocity  shock  wave  (Figure  3.1).  The  point  from  which  the 
jet  emanates  Is  a  stagnation  point  In  the  frame  of  reference  of  the 
shock  wave.  The  quasi-steady  analysis  shows  that  the  stagnation 
enthalpy  Is  approximately 


where  uA  Is  the  shock  velocity  In  the  frame  of  reference  of  the 
undisturbed  medium.  The  analysis  holds  only  for  a  constant  shock 
velocity  (steady  flow  In  the  shock  frame),  but  Equation  3.1  should  be 
a  reasonable  approximation  when  the  shock  speed  decays  gradually,  as 
In  the  spherical -source  experiments  of  PI. 
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Condensate  jetting  should  occur  as  long  as  h$  Is  sufficient  to 
melt  some  of  the  wall  Material:  1.e.t  h$£2KJ/g,  or  uA>2km/sec, 
for  steel  liners.  This  is  not  to  say  that  the  jet  will  remain  liquid. 
The  jetting  material  cools  as  It  flows  into  the  jet  and  is  likely  to 
become  small  solid  particles.  However,  when  h$  is  Insufficient  to 
cause  melting,  the  "Incompressible  fluid"  analysis  on  which  Equation 
3.1  Is  based  no  longer  holds  because  the  material  retains  considerable 
strength  In  the  plastic  state.  This  suggests  that  the  jetting  process 
will  probably  cease,  or  greatly  diminish,  when  Uj,  drops  below  about  2 
km/ sec. 


The  above  hypothesis  can  be  compared  with  experimental  results. 
Figure  3.2  shows  the  trajectories  of  the  ground  shock,  the  condensate 
jet,  and  the  vapor  arrival  (data  from  Ref.  2).  A  reference  line  at 
2  km/ sec  Is  also  shown  for  comparison.  The  ground  shock  velocity 
drops  below  this  value  at  an  axial  position  of  about  90  cm  (60  cm  from 
the  HE  source).  In  the  simplified  model  the  velocity  of  the  collapse 
point  (stagnation  point)  is  equated  with  the  ground  shock  velocity; 
for  the  actual  decaying  shock,  the  collapse  point  will  start  to  lag 
the  ground  shock  Increasingly  In  the  later  stages  of  collapse.  It  Is 
concluded  that  particulate  jetting  should  end  at  about  60  cm  from  the 
HE,  or  a  little  less. 

PI  experiments  have  provided  considerable  Insight  Into  the 
length  of  the  source  of  the  penetrating  Jet.  Use  of  "standoff" 
regions  of  wet  sand,  10  to  SO  cm  In  length,  between  the  HE  and  the 
beginning  of  the  pipe  resulted  In  steadily  diminishing  damage  to  the 
target.  As  shown  In  Figure  3.3,  the  extrapolation  of  the  results 
predicts  "zero  damage"  at  about  60  cm  standoff.  Other  experiments 
using  fu’l  length  pipes  with  Internal  helical  spirals  extending  from 
the  HE  only  for  part  of  the  pipe  length  (to  suppress  the  close-in  jet 
source,  as  In  the  standoff  cases)  produced  results  consistent  with  the 
standoff  results  (Fig.  3.3).  Thus  the  experimental  results  are  consis¬ 
tent  with  the  view  that  the  penetrating  jet  Is  formed  by  melting  of 
wall  material  on  collapse. 
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At  a  stagnation  enthalpy  somewhat  higher  than  the  melt  level  one 
expects  the  pipe  wall  to  vaporize  at  least  partially.  The  ground  shock 
velocity  of  4.2  km/ sec  near  the  source  should  produce  a  stagnation 
enthalpy  of  about  8.8  KJ/g  according  to  the  steady  state  theory.  Even 
larger  transient  values  are  seen  In  the  full  2D  calculation  at  the 
onset  of  pipe  collapse.  At  the  stagnation  point,  pressures  of  many 
hundreds  of  kbars  are  experienced  which  are  well  above  the  critical 
pressure  of  Iron  (thought  to  be  In  the  range  of  10  to  20  kbars).  The 
distinction  between  liquid  and  vapor  Is  meaningless  above  the  critical 
pressure.  However,  pressures  well  below  critical  occur  just  ahead  of 
the  stagnation  point  (according  to  the  calculations)  and  partial  vapor¬ 
ization  In  the  jet  Is  expected  for  enthalpies  of  3  to  4  KJ/g  or  more 
In  this  region.  Calculations  and  data  show  that  vapor  Is  a  significant 
component  of  the  jetting  material. 

The  above  arguments  suggest  that  the  Increased  target  damage 
observed  In  the  PI  Lead-Hrap  case  Is  due  In  part  to  a  reduction  of  the 
vapor  fraction  In  the  jetting  material.  The  calculated  total  energies 
for  the  jetting  vapor  plus  condensates  were  nearly  Identical  for  the 
two  cases. 

Since  the  penetrating  material  Is  liquid  or  solid,  but  plenty  of 
vapor  Is  also  formed,  the  question  arises  how  target  damage  might 
change  as  the  enthalpies  of  phase  change  are  varied.  At  first  glance, 
It  would  seem  desirable  to  seek  a  material  with  a  higher  heat  of 
fusion.  However,  this  normally  Implies  a  higher  heat  of  vaporization, 
and  some  energetic  vapor  generated  close  In  would  be  replaced  with 
condensates.  This  looks  like  a  way  of  exchanging  lower  velocity 
particulates  for  higher  velocity  particulates,  which  Is  the  wrong  way 
to  go.  If,  Instead,  one  seeks  a  lower  energy  of  vaporization,  the 
question  Is  whether  the  formation  of  more  vapor  comes  at  the  expense 
of  condensates,  rather  than  as  an  Increase  of  the  total  jetting.  An 
adverse  effect  Is  expected  to  occur:  material  that  Is  easily  vaporized 
tends  to  be  easily  melted  as  well,  and  the  condensate  jet  formation 
might  then  continue  to  considerably  lower  velocity  of  the  collapse 


point.  But  the  associated  specific  energy  ("penetrating  power"  per 
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unit  oass)  would  also  diminish  In  proportion  to  u^  .  In  effect, 
then,  the  use  of  wall  materials  with  lower  heats  of  phase  change  could 
reduce  target  damage  unless  the  reduction  of  specific  particulate 
energy  is  offset  by  a  proportionate  Increase  In  mass  of  jetting 
particulates.  This  question  deserves  careful  Investigation  from  both 
the  calculatlonal  and  experimental  point  of  view. 

The  practical  Implication  of  the  above  hypothesis  Is  that  it  may 
be  just  as  useful  to  reduce  the  fraction  of  jetting  energy  in  condensed 
states  as  to  reduce  the  total  jetting  energy.  Even  more  intriguing  Is 
the  fact  that  the  postulated  mechanism  for  the  successful  operation  of 
Internal  spiral  helices  -  turning  of  the  helix  Into  the  path  of  the 
particulate  jet  by  the  vapor  flow  -  would  imply  the  same  thing.  If  the 
jet  Is  Initially  all  vapor  and  the  condensate  jet  begins  later  and  at 
lower  speed,  there  Is  more  time  for  the  mechanism  to  work.  On  tapered 
pipes  this  could  become  crucial  (LS~5  showed  that  helical  spirals 
become  less  effective  as  the  taper  Increases). 

The  scenario  for  helix  operation,  while  still  tentative,  has 
been  made  more  credible  by  the  results  of  the  LS-5  experiment  which 
was  performed  after  the  publication  of  Ref.  1.  It  was  found  that 
loose  spirals  are  more  effective  than  attached  spirals.  It  was  also 
found  that  the  helix  Is  more  effective  than  a  series  of  rings.  In 
terms  of  the  aerodynamic  turning  scenario.  It  would  be  expected  that 
the  turning,  once  begun,  would  be  propagated  along  the  pipe  more 
readily  than  new  turning  could  be  Initiated  at  each  ring. 

In  summary*  a  phenomenological  description  of  the  mechanisms  of 
target  damage  by  jetting  and  of  the  prevention  of  target  damage  by 
Internal  spiral  helices  has  evolved.  These  are  consistent  with  the 
growing  body  of  LS  data.  The  practical  Implications  of  these 
suggested  mechanisms  are  extensive.  It  Is  Important  to  obtain  as  much 
Information  as  possible  on  LS-6  to  permit  a  sufficient  level  of 
confidence  In  either  these  or  alternative  scenarios. 
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3.3  TAPERED  PIPES 

The  LS-5  experiment  by  PI  Included  tapered  pipes  with  "nominal" 
LOS  taper  (0.375*  half  angle),  twice  nominal  taper,  and  four  times 
nominal  taper.  There  were  two  pipes  of  each  taper  without  a  helix; 
various  helical  and  Internal  ring  arrangements  were  also  used.  The 
following  table  summarizes  the  results  without  helices  and  compares 
them  with  the  standard  model  results. 


TABLE  3.1.  SUMMARY  OF  PI  DATA,  STANDARD  YS  TAPERED  PIPES 

(Refs.  2a,  b,  c,  d) 


Standard  (zero  taper) 

Nominal  Taper 
(0.375*  half  angle) 

Twice  Nominal  Taper 


No.  of  Samples  Mean  Hole  Volume  (cc) 


17.8  ±  8.3 


14.35±?  (11.0  and  17.7) 


Four  Times  Nominal  Taper 


17.55±?  (14.8  and  20.3) 
18. 8±?  (both  18.8) 


A  pre-shot  calculation  was  carried  out  with  the  STREAK  code  for 
the  4  x  nominal  taper  for  comparison  with  the  standard  model  calcula¬ 
tion.  Figure  3.4  shows  the  growth  of  energy  In  the  particulate  jet 
(Ej)  with  time  for  each;  Figure  3.5  shows  the  total  energy  In 
particulates  and  vapor.  Using  the  assumption  that  the  hole  volume  In 
the  target  Is  proportional  to  Ej,  as  In  Ref.  1,  one's  Initial 
Inclination  would  be  to  predict  that  the  hole  volume  In  the  target 
would  be  more  than  three  times  as  large  In  the  4  x  nominal  taper  model 
as  In  the  standard  model.  The  data  suggest  that  such  a  result  Is  not 
credible  for  any  realistic  values  of  the  uncertainty  In  hole  volume. 


Time  (usee) 


Figure  3.4.  Growth  of  kinetic  energy  in  the  particulate  jet 


Figure  3.5.  Growth  of  Total  Energy  (Particulate  jet  +  vapor) 


Further  consideration  of  the  calculatlonal  results  provides 
useful  Insight.  Notice  that  each  calculation  predicts  a  nearly-1 inear 
growth  of  Ej  and  E^^  for  a  time,  followed  by  a  gradual  decay  of 
growth  rate.  The  two  calculations  behave  differently  at  late  times. 
Energy  continues  to  Increase  In  the  taper  case  but  levels  off  In  the 
standard.  If  the  difference  In  the  late  stages  can  be  attributed  to 
numerics  as  discussed  below,  and  the  comparison  Is  regarded  as  mean¬ 
ingful  only  In  the  rapid  growth  stage,  the  calculatlonal  comparison 
becomes  much  more  consistent  with  the  data. 


Numerics  can  be  critical  at  the  stage  where  the  particulate 
jetting  Is  ending.  These  calculations  using  0.9  mm  radial  zoning  are, 
admittedly,  marginally  zoned  for  describing  the  jetting  of  collapsing 
0.3  nm  thick  pipe  walls  (the  fully-collapsed  pipe  extends  2.6  zones 
radially).  Such  zoning  was  necessary  for  economic  reasons  associated 
with  the  large  ratio  of  pipe  length  to  diameter.  The  comparison  of 
the  standard  and  Pb-Wrap  calculations  In  Ref.  l  Is  meaningful  only  In 
a  relative  sense  and  only  because  precisely  the  same  zoning  was  used 
in  each.  Ej  stopped  Increasing  In  each  after  only  about  20  cm  of 
pipe  was  collapsed,  whereas  the  actual  source  region  of  the  jetting 
appears  to  extend  to  about  60  cm  from  the  HE  source.  An  Intrinsic 
problem  with  the  tapered  pipe  calculation  is  that  the  jetting  source 
and  collapsing  pipe  cannot  both  have  exactly  the  same  zonal  resolution 
as  in  the  standard  model .  The  standard  used  10  zones  to  describe  the 
pipe  radius.  The  taper  calculation  accomplished  the  same  thing  by 
using  tapered  zoning  (this  Is  a  feature  of  STREAK'S  generalized 
coordinates);  In  this  case,  however,  the  collapsed  pipe  wall  Is 
resolved  In  fewer  zones  (Its  radius  fully  collapsed  Is  proportional  to 
the  square  root  of  Its  Initial  radius).  If  the  taper  of  the  zoning 
were  selected  to  preserve  the  same  resolution  of  the  collapsed  wall, 
then  the  pipe  radius  would  be  better  resolved  than  In  the  standard. 
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Short  of  zoning  both  calculations  more  finely,  all  that  can  be 
done  Is  to  compare  the  rates  of  growth  of  In  the  region  before 
the  particulate  jet  source  becomes  too  weak  to  be  sustained  with  the 
zonal  resolution  employed.  The  calculated  ratio  of  tapered/standard 
Ej's  Is  then  Interpreted,  roughly,  as  1.6±0.3. 

The  above  Interpretation  Is  not  at  all  Inconsistent  with  the 
data.  There  Is  an  Indication  In  the  data  of  a  trend  toward  Increasing 
hole  volume  with  Increasing  taper.  Unfortunately,  the  results  at  zero 
taper  do  not  fit  this  trend.  However,  note  that  the  four  standard 
models  of  LS-3  seem  anomalously  high  (see  Section  3.4).  If  the  LS-3 
results  are  due  to  some  unknown  Inconsistency  between  the  driving 
conditions  In  that  experiment  and  the  others,  and  therefore  could  be 
Ignored,  V  would  become  12.35  cc  for  the  standard  which  fits  the 
observed  trend  perfectly.  This  exercise  would  produce  a  ratio  of 
V  (4  x  nominal  taper)/V  (standard) =  1.5  which  Is  totally  consistent 
with  the  calculations. 

Unfortunately,  there  Is  no  substantial  reason  for  Ignoring  the 
LS-3  standard  results  at  this  time.  All  that  can  be  said  Is  that  the 
calculatlonal  results,  prior  to  the  rapid  decay  In  energy  growth,  are 
not  denied  by  the  data  because  of  significant  uncertainty  In  V  for  the 
tapered  pipes.  It  should  also  be  pointed  out  that  this  revised 
Interpretation  of  the  calculations  does  not  change  the  comparison  of 
lead  wrap  vs  standard  (Ref.  1 )  because  the  slopes  In  the  rapid  growth 
stage  In  those  calculations  are  proportional  to  the  peak  values  of 


One  further  comment  about  the  calculatlonal  results.  The 
calculations  predict  approximately  equal  total  energy  (vapor  plus 
particulates)  for  the  Standard  and  taper  cases  but  the  taper  case  has 
less  energy  In  vapor  and  more  In  particulates  than  the  standard.  This 
Is  qualitatively  consistent  with  available  tlme-of-arrlval  data  for 
the  vapor  (Ionization  pins).  The  standard  models  (LS-2,3,4)  Indicated 


vapor  arrival  velocities  of  about  0.85  to  1.2  cm/us  whereas  the  only 
tapered  pipe  (without  helix)  reporting  data,  twice  normal  taper. 
Indicated  0.79  cm/us  (Ref.  2d).  This  Is  Indicative  that  the  vapor 
flow  Is  less  energetic  In  a  tapered  pipe. 

3.4  STANDARD  MODEL  STATISTICS,  REVISITED 

Additional  Insights  have  been  gained  by  reviewing  the  set  of 
eight  Standard  model  results  for  the  hole  volume  In  the  target. 

First,  It  has  been  noticed  that  a  lognormal  distribution  fits  the  data 
better  than  does  a  normal  distribution.  This  Is  seen  In  Figure  3.6. 
This  observation,  previously  overlooked.  Is  completely  consistent  with 
the  known  tendency  of  quantities  which  cannot  take  negative  values. 
(The  normal  distribution  predicts  that  about  2%  of  Standard  models 
will  have  negative  hole  volumes.) 

It  has  also  been  noticed  that  the  data  spread  Is  narrower  If  the 

eight  samples  of  Table  3.2  are  split  Into  two  distinct  sets:  the 

four  samples  from  LS-3  and  the  other  four  from  LS-2,4.  The  actual 

relative  sample  variance  (S)  Is  reduced  even  more  by  the  separation 

Into  two  samples  (the  "unbiased"  standard  deviation,  o,  given  by 
2  2 

o  *  S  (N/N-l)  for  a  sample  of  N  results.  Is  more  Inflated  for 
smaller  samples).  The  values  of  S/V  are  0.31  for  LS-3,  0.27  for  the 
other  four,  and  0.44  for  the  combined  eight.  This  suggests  a  likeli¬ 
hood  that  LS-3  was  different  from  the  others  In  some  way.  The  shock 
trajectory  In  the  wet  sand  was  several  percent  faster  In  LS-3  than  In 
LS-2  (Ref.  2c),  which  In  Itself  doesn't  seem  a  large  enough  effect  to 
explain  a  factor  of  two  In  V,  but  some  related  difference  In  the 
preparation  of  the  medium  or  In  the  HE  detonation  may  explain  It.  The 
separation  of  the  standard  experiments  of  LS-3  from  those  of  LS-2,4 
remains  pure  speculation,  and  even  If  It  could  be  justified.  It  would 
not  provide  the  much-desired  justification  for  neglecting  the  LS-3 
results  In  the  analysis  of  the  LS-5  tapered  pipe  results  (Section 
3.3):  It  seems  the  ground  shock  was  just  as  fast  In  LS-5  as  In  LS-3 
(Ref.  2d). 


SECTION  4 
LOS  FLOW 

This  section  discusses  recent  Improvements  In  the  phenomeno¬ 
logical  understanding  of  LOS  flows  and  reviews  the  state-of-the-art  In 
their  calculatlonal  modeling  (Section  4.1).  It  also  discusses  the 
Implications  of  recent  results  for  the  optimization  of  closure  design 
(Section  4.2).  Some  Investigations  carried  out  for  specific  shots  are 
described  In  Section  4.3. 

4.1  PHENOMENOLOGY  AND  MODELING 

The  two  most  recent  horizontal  LOS  events,  Diablo  Hawk  (DH)  and 
Miners  Iron  (MI),  produced  LOS  diagnostic  data  of  a  quality  not 
previously  achieved.  Each  of  these  events  produced  an  Internally 
consistent  set  of  data  at  a  number  of  stations  along  th^  LOS  and  the 
results  of  the  two  events  were  consistent  with  each  other.  This  data 

base  makes  It  possible,  for  the  first  time,  to  test  hypotheses  about  _ 

the  phenomenology  of  LOS  flow. 

The  pressure  history  data  from  both  events  strongly  Indicate  two 

distinct  LOS  pressure  pulses  at  each  station  prior  to  ground  shock 

arrival.  This  Is  seen  In  the  examples  given  in  Figures  4.1  and  4.2; 

data  from  all  other  stations  show  similar  behavior.  Figures  4.3  and 

4.4  show  peak  pressure  vs  range  for  the  first  and  second  pulses, 

respectively,  from  DH  and  MI.*  A  striking  feature  of  the  results  is 

the  qualitative  difference  In  the  decay  of  pressure  with  range 

-3 

associated  with  the  two  pulses.  The  first  pulse  decays  as  R 
whereas  the  second  manifests  an  approximately  exponential  decay  with 
range.  (The  upturn  In  the  observed  amplitude  trends  just  before  the 
MACs  In  Figure  4.4  Is  probably  caused  by  the  encounter  with  the 
detonation  products  and  nitrogen  released  during  the  MAC  actuation.) 


*  The  data  were  provided  by  R.  C.  Bass  and  C.  W.  Cook  of  Sandla 
National  Laboratories. 
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These  data  challenge  the  analyst  to  provide  phenomenological 
explanations.  The  rest  of  this  Section  Is  devoted  to  a  presentation  of 
such  explanations  and  the  reconciliation  of  calculations  with  the  data. 

Previous  experience  with  LOS  closure  calculations  provides  ample 
evidence  of  two  pulses.  The  first  pulse  Is  attributed  to  prompt  energy 
Injection  Into  the  LOS,  primarily  In  the  blowoff  gases;  this  Is  formed 
by  radiation-induced  vaporization  of  the  pipe  wall  and  has  no  analog  In 
the  HE-driven  simulations  by  PI.  The  second  pulse  Is  a  consequence  of 
further  energy  Injection  Into  the  LOS  by  action  of  the  ground  shock; 
this  is  analogous  to  the  vapor  flow  In  the  PI  simulations  with  one 
important  difference:  in  the  simulations  this  pulse  Is  an  expansion 
into  a  vacuum  whereas  In  the  LOS  it  Is  a  shock  propagating  Into  the 
residual  gases  of  the  first  pulse.  This  probably  explains  the  observed 
difference  in  velocity.  The  calculatlonal  comparisons  presented  later 
In  this  Section  support  this  view. 

The  Initial  pressure  associated  with  the  blowoff  gases  decreases 

with  range,  and  the  associated  pressure  gradient  rapidly  Initiates  a 

flow  in  the  pipe.  The  more  dense  blowoff  gases  nearer  the  source  begin 

to  expand  Into  the  near-vacuum  farther  from  the  source.  The  expansion 

Is  a  combination  of  radial  expansion  resulting  from  the  conical  shape  of 

the  pipe  and  axial  expansion  along  the  pipe  direction.  Without 

significant  wall  Interaction,  this  expansion  would  result  In  the  one-way 

conversion  of  plasma  Internal  energy  into  directed  kinetic  energy  and  an 

-y 

associated  pressure  drop.  The  pressure  would  decay  as  Vpiasma,  where 
V  Is  the  volume  of  a  Lagranglan  element  of  the  blowoff  plasma. 

However,  LOS  flow  calculations  have  always  indicated  much  too 
early  arrivals  of  pressure  pulses  along  the  axis  unless  ablation  Is 
considered.  Ablation  of  the  wall  continually  adds  mass  to  the  flow 
and  dissipates  kinetic  energy.  Calculations  considering  only  the 
blowoff  gases  (discussed  below)  Indicate  that  the  partition  between 
directed  and  Internal  energy  remains  nearly  constant  after  an  initial 
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transient  behavior.  That  Is,  the  dissipation  and  expansion  effects 
offset  each  other.  In  this  case,  the  pressure  of  an  element  should 
vary  Inversely  as  its  volume: 

poc  y_1 

As  the  element  moves  to  Increasing  range  (R)  its  cross-section  varies 
as  R  because  of  the  conical  pipe  and  Its  axial  length  may  be 
assumed  to  vary  as  R~^1+e^  where  e  may  be  positive  or  negative  and 
Is  thought  to  be  small.  This  Implies  a  decay  of  pressure  with  range 
of  the  form 

dP/P  °c  -  3dR/R , 

P/?n  =  (R/R  )"3  (4.1) 

o  o 

which  is  consistant  with  the  experimental  data  except  at  close  range. 

A  refinement  of  the  above  approach  takes  Into  account  the  fact 
that  the  pipe  diameter  does  not  approach  zero  with  R.  One  alternative 
formulation  Is 

dP/P  »  -  3dD/D 

where  the  pipe  diameter,  D,  is  given  by 

D  -  D  +  <*(R-R  )  . 
o  o 

Ro*  D0  are  va*ues  at  some  r®ference  point.  Integration 
produces  the  following: 


Another  alternative,  seemingly  more  defensible,  superposes  an  axial 
expansion  term,  -dR/R,  and  a  radial  expansion  term,  -2dD/D,  to  obtain 
dP/P.  The  result  of  integration  Is: 

P/P0  =  <R/Ro)-1  (D/D0)"2  (4.3) 

The  predictions  of  Eqs.  (4.2)  and  (4.3)  for  decay  of  the  first  pulse 
are  shown  in  Figure  4.3.  (Each  curve  was  fitted  to  the  data  at  one 
point.)  It  Is  seen  that  both  are  superior  to  Eq.  (4.1)  for  small  R, 
with  Eq.  (4.3)  slightly  favored,  whereas  all  three  are  equivalent  at 
long  range  (where  dD/D  =  dR/R).  The  parameters  used  In  producing  the 
curves  are  shown  In  Table  4.1 


TABLE  4.1.  LOS  PIPE  PARAMETERS  FOR  DIABLO  HAWK  AND  MINERS  IRON 


DH  MI 


a  0.01309  0.0144 

R  (m  from  WP)  11  11 

o 

D„  (m),  at  R  0.305  0.335 

o  o 


The  second  pulse  is  regarded  as  a  shock  wave  because  it 
propagates  Into  the  residual  gas  of  the  first  pulse.  Extensive 
experimental  Information  on  the  decay  of  strong  shock  waves  in  straight 
pipes  exists,  e.g.,  the  material  sunmarized  In  Ref.  (5)  (both  nuclear 
and  HE  driven).  An  exponential  deciy  with  range  characterizes  such 
shocks  when  the  flow  behind  the  shock  Is  hot  enough  to  ablate  the  wall: 

P  «  exp(-KR/D) 

•P 
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where  K  is  a  constant.  K=  0.0576  Is  consistent  with  most  of  the  data. 
In  particular,  from  the  nuclear  shock  tube  events,  Hybla  Gold  and 
Marvel.  The  R/D  dependence  can  be  justified  analytically.  The 
observed  end  of  the  exponential  decay  regime  Is  also  consistent  with 
the  postulated  ablation  process.  For  shocks  Into  ambient  air,  the 
exponential  decay  gives  way  to  a  more  gradual  decay  at  just  the  point 
where  the  post-shock  air  Is  no  longer  hot  enough  to  vaporize  the  wall 
material  (P$~  100  to  300  bars).  For  shocks  Into  a  residual  plasma, 
say  Iron  vapor,  ablation  should  play  a  role  down  to  considerably  lower 
pressures. 

In  a  tapered  pipe.  It  Is  proposed  that  the  decay  of  shock 
pressure  with  range  be  expressed  as  a  superposition  of  ablation- 
dominated  shock  decay  as  described  above  and  of  taper  effects  as  In 
the  analysis  of  the  first  pulse.  The  analog  of  Eq.  (4.2)  Is  obtained 
by  setting 

dP/P  -  -  KdR/D  -  3dD/D 

Substitution  of  the  pipe  taper  relation,  followed  by  Integration, 
produces 

P/PQ  *  (D/D^)”^ 3+K/a)  (4.4) 

Alternatively,  the  analog  of  Eq.  (4.3)  Is  obtained  In  similar  fashion: 

P/PQ  -  (R/Rq)-1  (D/D0)"(2+K/a)  ^4'5) 

The  predictions  of  Eqs.  (4.4)  and  (4.5)  for  the  decay  of  the  second 
pulse  are  shown  In  Figure  4.4.  The  predictions  are  generally  quite 
consistent  with  the  data.  Eq.  (4.4)  appears  better-behaved  for  small 
R;  otherwise,  It  Is  difficult  to  choose  between  the  two.  (The  analog 
of  Eq.  (4.1)  behaves  worse  at  small  R  and  Is  not  shown.) 
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The  above  exercise  Indicates  that  the  physical  assumptions  on 
which  It  Is  based  are  reasonable:  namely,  that  the  first  pulse  Is  an 
expansion  Into  a  near-vacuum,  the  second  pulse  Is  a  shock  wave,  and  both 
are  strongly  Influenced  by  ablation.  There  Is  some  uncertainty  as  to 
details,  but  the  qualitative  behavior  of  all  of  the  alternative  expres¬ 
sions  Is  reasonable.  It  Is  especially  notewothy  that  the  quasl-analytlc 
expressions  derived  here  are  consistent  with  the  experimental  observation 
that  the  amplitude  of  the  second  pulse  exceeds  the  first  at  short  range 
but  the  opposite  Is  true  at  long  range. 

The  Diamond  Ace  event  has  been  designed  with  a  larger  pipe  taper 
than  usual.  Therefore,  It  constitutes  a  useful  test  of  the  validity  of 
the  above  quasl-analytlc  techniques  of  as  well  as  of  calculatlonal  tech¬ 
niques  for  LOS  flow.  Investigations  of  the  effect  of  LOS  taper  on  second 
pulse  decay,  using  Equations  4.4  and  4.5,  have  shown  that  although  the 
term  D  In  the  argument  grows  more  rapidly  If  a  Is  Increased,  this  Is 
somewhat  offset  by  the  accompanying  reduction  of  the  exponent.  That  Is, 
the  ablation-dominated  aspect  of  shock  decay  becomes  less  prominent  with 
Increased  taper.  Nevertheless,  an  Increase  of  taper  for  fixed  values  of 

R,  R  ,  0  and  K  does  result  In  a  reduction  of  P/P  for  all  combina- 
oo  o 

tlons  likely  to  be  Interesting.  We  look  forward  to  the  comparison  of 
Diamond  Ace  data  with  computational  and  quasl-analytlc  results. 

A  brief  comment  on  muffler  effects  Is  In  order.  The  muffler 
location  Is  Indicated  In  Figures  4.3  and  4.4  (approximately  the  same  for 
both  events).  There  Is  no  Indication  of  a  change  In  the  trend  of 
pressure  decay  at  the  muffler.  We  Interpret  this  as  evidence  that  the 
muffler  has  little  or  no  effect  on  flow  amplitudes.  However,  R.  C.  Bass 
of  SNLA  Interprets  the  data  somewhat  differently  and  believes  the 
muffler  has  some  effect. 

A  comparison  of  calculatlonal  predictions  with  the  MI/DH  data  set 
Is  In  order.  The  desired  calculations  would  be  STREAK  calculations 
Initiated  from  RAD9  calculations  as  early  and  as  correctly  as  possible. 
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This  Is  not  quite  what  exists,  as  we  shall  see,  but  what  exists  Is  ade¬ 
quate  for  an  assessment  of  how  well  the  current  calculatlonal  technique 
performs. 

A  preshot  LOS  flow  calculation  for  Diablo  Hawk  was  performed  in 
1976  with  the  UNION  code.  The  FLIP  segment  then  In  use  employed  an 
ablation  model  that  Is  now  obsolete.  The  calculation  was  Initiated  at 
0.1  ms  from  a  RAD9  solution  and  carried  to  4  ms.  The  primary  purpose  of 
that  calculation  (which  was  the  last  In  the  series  performed  as  part  of 
the  DH  closure  design  effort)  was  a  comparison  of  predictions  of  Injec¬ 
ted  plasma  and  particulate  energy  Into  the  LOS  with  those  predicted  for 
other  configurations.  The  secondary  purpose  was  to  make  preshot  LOS 
pressure  history  predictions  to  compare  with  the  eventual  data;  for 
this  purpose  an  LOS  calculation  was  continued  on  stand-alone  FLIP  from  4 
ms  to  more  than  60  ms.  The  calculated  arrivals  proved  to  be  quite  late 
compared  with  the  measurements. 

Later,  the  ablation  modeling  was  completely  overhauled  following 
Hybla  Gold.  The  revised  model  was  shown  to  give  much  better  agreement 
with  observed  arrivals  -  at  least  to  the  muffler  (see  a  classified 
report  on  LOS  studies,  1981);  the  less  rapid  ablation  obtained  with  the 
new  model  allowed  faster  arrivals.  STREAK  then  replaced  UNION  as  the 
primary  20  code  for  closure  calculations.  It  employed  the  up-to-date 
ablation  modeling  for  pipe  flow.  The  DH  closure  calculation  was  repeated 
from  0.1  to  4  ms  with  STREAK  In  1981.  The  LOS  flow  calculation  after 
4  ms  was  again  continued  with  the  stand-alone  FLIP.  Although  the  calcu¬ 
lated  plasma  arrivals  were  much  Improved,  they  were  still  distinctly 
later  than  observed. 

The  results  of  the  STREAK/FLIP  calculation  and  the  experimental 
data  shown  In  Figure  4.5  for  the  42.3  m  station  are  typical.  The 
calculated  first  pulse  begins  at  about  4  ms  whereas  the  observed  first 
pulse  begins  at  about  2.5  ms.  The  calculated  second  pulse  begins  at 
about  8.5  to  9  ms,  less  than  1  ms  late  -  a  much  smaller  percentage  error 
in  pulse  speed  than  with  the  first  pulse.  Both  calculated  pulse 
amplitudes  are  too  large  by  more  than  a  factor  of  two. 
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7 THE.  IN  SECGN03  FROM  ZERO  FTC'J 

- 2D  STREAK  to  4  ms,  continued  on  FLIP  (Quasi-ID,  see  Refs.  2,8) 

Blowoff  close-in  onlyfsee  text) 

-  Blowoff  only,  FLIP 

- -  Expected  result  if  full  blowoff  included  in  the  full  2D  calculation. 


Figure  4.5.  Calculated  and  measured  LOS  pressure  histories 
from  Diablo  Hawk. 
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The  reason  for  the  discrepancy  In  the  first  pulse  onset  time  was 
found  to  be  In  the  treatment  of  the  prompt  blowoff  vapor.  The  RAD9 
calculation  of  the  first  0.1  ms  had  Included  prompt  blowoff  only  In  the 
Initial  grid  which  extended  a  short  distance  beyond  the  zero  room. 

As  that  calculation  progressed  and  the  spatial  grid  was  expanded,  the 
added  sections  of  LOS  were  Initiated  with  nominal  low-density,  low- 
pressure  vapor  conditions.  This  was  of  no  real  significance  to  the 
determination  of  LOS  closure  behavior  or  of  energy  Injected  Into  the 
LOS.  However,  It  resulted  In  the  loss  from  the  calculation  of  signifi¬ 
cant  Initial  blowoff  vapor  as  well  as  additional  LOS  vapor  which  left 
the  grid  during  each  stage  of  calculatlonal  zoning. 

We  decided  to  Initiate  a  FLIP  calculation  with  a  complete  repre¬ 
sentation  of  the  OH  blowoff  vapor.  Although  no  additional  energy 
Injection  was  considered,  this  was  expected  to  provide  a  fair  represen¬ 
tation  of  first  pulse  arrival.  The  comparison  of  the  results  with  the 
first  pulse  onset  data.  Figure  4.6,  supports  this  view  completely.  The 
wave  form  at  42.3  m  from  this  calculation  Is  shown  In  Figure  4.5;  also 
shown  Is  an  estimate  of  how  this  first  pulse  would  peak  If  It  were 
determined  In  a  full  2D  calculation  describing  the  continuing  LOS  energy 
Injection  process.  The  results  Indicate  that  the  present  calculatlonal 
tools  are  highly  credible  as  to  their  ability  to  predict  two  distinct 
LOS  plasma  pulses  with  appropriate  arrival  times.  Further  confirmation 
of  this  claim  Is  given  In  the  comparison  of  calculated  and  measured 
second  pulse  arrivals  for  DH  {Figure  4.6).  A  blowoff-only  calculation 
has  also  been  performed  for  Miners  Iron.  The  calculated  arrival  times 
are  again  In  good  agreement  with  observed  first  pulse  arrivals,  (Figure 
4.6). 


The  calculated  pressures  In  Figure  4.5  are  2  to  3  times  too  high. 
The  same  was  found  at  other  stations  In  Diablo  Hawk.  This  chronic 
problem  with  LOS  calculations  Is  attributed  to  the  Inability  of  the 
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Figure  4.6.  Comparison  of  measured  and  calculated  pulse  arrivals.  Miners  Iron 
and  Diablo  Hawk. 


current  modeling  to  treat  quenching  mechanisms  such  as  massive  entrain¬ 
ment  of  exposed  grout  after  pipe  rupture  during  expansion  and  loss  of 
vapor  Into  cracks  and  Into  the  advancing  ground  shock-driven  collapse 
region. 

Despite  this  shortcoming.  It  Is  apparent  that  calculatlonal 
techniques  have  Improved  substantially  during  the  past  several  years. 

The  ablation  modeling  enjoys  growing  credibility  by  virtue  of  calcula¬ 
tlonal  results  that  are  consistent  with  a  broad  data  base  (see  also  the 
comparisons  presented  In  Ref.  (5)).  Admittedly,  the  evidence  for  abla¬ 
tion  Is  only  Indirect.  However,  other  calculations  have  shown  that  when 
ablation  Is  Ignored  In  high-enthalpy  flows, the  calculated  pressure 
pulses  arrive  much  too  early  and  have  much  too  low  amplitude.  And  it  is 
especially  compelling  that  the  rate  of  attenuation  of  a  shock  propagating 
Into  an  air-filled  pipe  has  been  observed  to  undergo  a  significant  reduc¬ 
tion  at  just  the  range  where  the  air  temperature  behind  the  shock  becomes 
too  low  to  vaporize  the  wall  (Hybla  Gold  1  ft  pipe.  Ref.  5). 

It  Is  concluded  that  our  calculatlonal  tools  are  steadily 
Improving  with  the  aid  of  a  growing  base  of  high  quality  data  and  a 
persistent  effort  to  Improve  the  physical  modeling.  These  tools  are 
quite  credible  as  to  their  ability  to  predict  two  distinct  LOS  plasma 
pulses  with  reasonably  accurate  arrival  times.  The  tendency  to  over¬ 
estimate  pressures  Is  at  least  conservative  and  the  reasons  for  It 
appear  to  be  understood.  To  complete  the  analogy  with  the  PI  simula¬ 
tions,  the  2D  closure  calculations  for  the  LOS  also  indicate  the  forma¬ 
tion  of  a  jet  of  particulates.  The  leading  edge  of  this  Is  somewhat 
slower  than  the  calculated  second  shock.  There  are  insufficient  experi¬ 
mental  data  for  quantitative  comparison  with  calculations,  however, 
cratering  of  MAC/DAC  doors  has  been  observed  in  all  cases  where  those 
doors  have  been  recovered. 

4.2  CLOSURE  OPTIMIZATION 

It  Is  useful  to  review  the  position  of  closure  design  on  Its 
evolutionary  path  and  consider  what  Improvements  might  be  realizable  In 

the  near  future.  We  will  consider:  tuballoy  extensions;  helical  wraps; 
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a  proposed  high  explosive-driven  closure  FAC  (fast-acting  closure)  simi¬ 
lar  to  the  super  HE  machine  used  on  VLOS  events;  and  stemming  materials. 

The  development  and  evolution  of  tubal loy  extensions  over  four  DNA 
events  has  been  the  most  recent  major  change  in  closure  design.  This 
evolved  from  a  1975  calculatlonal  study.  There  is  strong  evidence  that 
their  use  has  reduced  the  Injection  of  energetic  vapor  into  the  LOS. 

This  Is  demonstrated  In  Figure  4.7  which  shows  measured  peak  pressure 
(the  larger  of  the  first  and  second  pulse)  vs  scaled  range  for  a  number 
of  events,  seven  with  Iron  extensions  and  three  with  tubal loy.*  Even 
Miners  Iron,  which  surprisingly  produced  peak  pressures  several  times 
higher  than  Diablo  Hawk  (see  also  Figures  4.3  and  4.4),  is  seen  to  be 
quite  benign  In  this  context.  The  Miners  Iron  result  Is  still  under 
study,  but  the  cause  of  the  unexpectedly  high  pressures  Is  thought  to  be 
the  larger  LOS  diameter  near  the  source,  not  the  performance  of  the 
tubal loy  extension. 

The  tubal loy  extension  reduces  energy  Injection  Into  the  LOS 
because  Its  high  density  slows  the  velocity  of  convergence  of  material 
toward  the  pipe  axis.  The  role  of  a  slower  extension  collapse  In  genera¬ 
ting  cratering  flow  must  be  carefully  considered.  As  discussed  In  Section 
3.2,  the  collapse  of  a  steel  liner  at  velocities  between  about  2  and  4 
km/s  should  produce  a  predominantly  condensate  jet.  The  Increased  target 
damage  In  the  PI  lead  wrap  simulation  may  be  due  In  part  to  the  produc¬ 
tion  of  a  jet  that  Is  less-fully  vaporized.  The  optimization  of  the 
tuballoy  extension  must  Include  careful  analysis  of  this  problem  of 
vapor/condensate  partition.  The  experience  with  PI  simulations  teaches 
that  such  analysis  requires  calculations  which  adequately  resolve  impor¬ 
tant  features  such  as  the  steel  liner  Inside  the  tuballoy.  Such  calcula¬ 
tions  were  begun  during  the  contract  period  but  did  not  reach  conclusive 
stages.  These  are  continuing  as  part  of  the  follow-on  research  effort. 


*  The  data  summarized  In  Figure  4.7  were  provided  mostly  by  R.C. 
Bass  of  Sandla,  though  some  other  sources  were  used  and  there  has  been 
some  Interpretation  by  the  authors. 
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The  next  Improvement  planned  for  closure  design  Is  the  helical 
wrap.  Analysis  of  laboratory  simulation  of  Its  performance  (see  discus¬ 
sion  In  Ref.  (1)  and  summary  In  Section  3)  Indicates  that  It  reduces 
target  damage  by  deflecting  or  intercepting  the  jetting  particulates. 

It  Is  currently  thought  (by  us  at  least)  that  the  helix  performs  this 
function  by  being  turned -into  the  path  of  the  particulates  by  the  faster- 
arriving  vapor  flow.  Because  of  this  apparent  need  for  some  vapor  flow 
and  the  strong  possibility  that  cratering  (particulate)  flow  Is  a  greater 
threat  to  experiments  and  containment  hardware  than  the  vapor  flows 
currently  experienced,  it  may  not  be  fruitful  to  try  to  effect  drastic 
further  reduction  of  vapor  Injection  by  extension  design.  It  may  even  be 
necessary  to  consider  making  the  helix  of  lower  density  material  when 
large  pipe  tapers  are  used,  to  enable  It  to  turn  into  the  stream  quickly 
enough  (note  results  obtained  In  Pi's  LS-5,  Ref  (2d). 

The  next  improvement  following  the  helical  wrap  may  be  an  economi¬ 
cal  modification  of  the  FAC  concept.  Some  preliminary  thinking  was  done 
during  the  contract  period  on  the  use  of  a  conical  device  that  Is  given 
a  modest  Initial  explosive- Induced  motion  toward  axis.  The  explosive 
could  be  detonated  by  prompt  radiation,  eliminating  detonation  Initiation 
devices  entirely.  LOS  plasma  would  be  relied  on  to  complete  the  closure. 
The  plasma  would  move  mostly  into  the  space  opened  around  the  cone  by  the 
HE  detonation  and  should  ultimately  make  the  closure  tight.  The  fact 
that  some  of  the  al ready-modest  LOS  plasma  energy  might  get  past  the 
closure  is  not  considered  significant.  Experience  with  PI  simulations. 
Indicating  that  jetting  particulates  move  preferentially  along  the  pipe 
walls,  suggests  that  such  a  device  could  Intercept  nearly  all  of  the  par¬ 
ticulates  even  If  not  fully  closed  when  they  arrive.  Such  a  device  offers 
the  possibility  of  many  economies.  For  example,  it  appears  capable  of 
doing  the  muffler's  job  even  better;  It  could  work  In  conjunction  with 
a  helical  wrap  but  could  also  replace  it;  and,  most  importantly.  It 
could  possibly  even  allow  at  least  one  of  the  MACs/DACs  to  be  eliminated. 
More  detailed  analysis  and  calculations  relating  to  this  concept  are 
being  carried  out  as  part  of  the  follow-on  research  effort.  The  central 
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problem  appears  to  be  the  need  to  Insure  that  the  plasma  will  not  cause 
such  energetic  collapse  of  this  configuration  that  particulate  jetting 
Is  Initiated.  This  would  merely  transfer  the  jetting  problem  from  its 
present  location  to  a  more  remote  location.  This  requirement  also 
Implies  that  such  a  closure  must  be  sufficiently  far  from  the  source 
that  ground  shock-driven  pipe  collapse  beyond  the  closure  is  not  a  source 
of  particulate  jetting.  This  may  require  that  It  be  located  about  30  m 
or  more  from  the  source. 

The  selection  of  stemming  materials  has  also  received  considera¬ 
tion.  The  Huron  King  experience  stimulated  reconsideration  of  the  use 
of  dry  sand  stemming  along  some  part  of  the  LOS  because  Its  high  air 
void  content  results  in  rapid  attenuation  of  the  ground  shock.  This  Is 
the  subject  of  a  STREAK  calculatlonal  study  begun  during  the  contract 
period.  The  two  configurations  compared  were  Miners  Iron  as-built  and 
with  sand  stemnlng  replacing  the  grout  beginning  1  m  beyond  the  end 
wall.  Several  cost-saving  approximations  were  used:  a  plugged  LOS  out 
to  the  end  wall;  simplification  of  the  material  representation  in  the 
zero  room;  and  omission  of  radiation  transport.  These  are  not  thought 
to  affect  the  questions  of  Interest  In  any  significant  way.  These 
calculations  are  part  of  the  study  mentioned  above  under  extension- 
optimization.  The  1  m  offset  was  selected  because  of  results  obtained 
with  an  earlier  RAD9  calculation.  RAD9  showed  an  Initially  more  rapid 
energy  Injection  into  the  LOS,  during  the  first  1  m  of  ground  shock 
propagation,  because  of  the  delay  In  closing  the  extension.  The  subse¬ 
quent  energy  Injection  rate  was  lower  than  with  grout  because  of  the 
weaker  shock  In  the  sand.  The  calculations  did  not  get  far  enough  In 
time  during  the  contract  period  to  allow  firm  conclusions  (120  ms); 
however,  the  results  Indicate  that  there  may  be  some  significant  (but 
not  dramatic)  reductions  In  LOS  energy  to  be  realized  from  optimization 
of  the  stemming  materials.  The  results  will  be  reported  In  detail  when 
the  study  Is  completed. 

Closure  optimization  may  be  entering  a  stage  of  rapid  development. 
There  Is  obviously  the  potential  for  great  Improvement  in  closure  perfor¬ 
mance  and  substantial  cost  reduction.  New  concepts  are  being  tried  and 


others  are  In  conceptual  development.  The  key  role  of  both  simulation 
experiments  and  calculatlonal  analysis  In  the  evolution  of  concepts  is 
obvious.  We  anticipate  rapid  progress  in  this  area  in  the  current 
research  effort. 
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SECTION  5 

RESIDUAL  STRESS  STUDIES 


5.1  SRI  GROUT  SPHERES 

SRI  International  has  since  1976  been  conducting  small  scale 
laboratory  experiments  to  study  the  residual  stress  fields  around 
explosively  formed  cavities  In  12  Inch  diameter  spheres  of  rock-matching 
grout.  S-CUBED  has  been  conducting  finite  difference  simulations  of  the 
high  explosive  detonation  and  the  subsequent  non-linear  dynamic  processes 
which  result  In  the  formation  of  compressive  residual  stress  fields  in 
the  grout  surrounding  the  explosive  cavity  in  order  to  validate  our 
capability  to  calculate  containment  related  phenomena.  Before  1981  the 
only  laboratory  data  available  to  be  directly  compared  to  the  dynamic 
calculations  consisted  of  final  cavity  radius  and  the  pressure  and 
integrated  impulse  at  a  quartz  gauge  emplaced  in  the  bottom  of  the  water 
tank  surrounding  the  grout  sphere.  In  1981,  SRI  began  to  emplace  particle 
velocity  gauges  inside  the  grout  spheres.  These  consisted  of  concentric 
current-carrying  loops  of  wire  in  a  magnetic  field  cast  symmetrically 
about  the  explosive  charge.  Here  we  discuss  the  constitutive  modeling 
advances  made  to  date  to  best  simulate  the  particle  velocity  records  in 
2C4  grout. 

Earlier  work  by  Rimer  and  Lie,  Reference  (6),  had  shown  that 
increases  in  strength  above  the  static  failure  surface  measured  in  the 
laboratory  were  necessary  to  match  the  measured  cavity  radius.  These 
increases  in  strength  were  assumed  to  be  proportional  to  the  logarithm 
of  strain  rate  experienced  by  a  grout  element.  Strain  rates  of  105  - 
107  s-1  were  calculated  at  the  shock  wave  for  these  small  scale 
experiments  so  that  the  increases  in  strength  were  quite  large.  The 
calculations  with  strain-rate-dependent  strength  gave  residual  stress 
peaks  of  almost  1.0  Kbar,  much  greater  than  the  0.40  -  0.45  Kbar 
fracture  initiation  pressures  measured  iy  hydrofracturing  from  the 
explosively  formed  cavity. 


Comparisons  made  between  the  measured  particle  velocity 
records  and  those  calculated  using  the  above  models  showed  two  major 
differences;  the  calculated  positive  velocity  pulses  were  somewhat 
narrower  than  the  measurements  and  the  negative  pulses  were  a  factor 
of  three  too  narrow.  An  extensive  effort  was  initiated  to  Improve  the 
constitutive  models,  and.  In  particular,  to  widen  the  negative  velocity 
pulse  (the  cavity  rebound  phase  of  the  phenomenology)  so  important  to 
the  calculation  of  compressive  residual  stress  fields.  The  constitu¬ 
tive  modeling  variations  examined  in  1981  are  discussed  beginning  with 
the  earliest  changes  and  ending  with  the  most  up  to  date. 

1.  Crush  curve  variations:  These  mainly  impacted  the  peak 
velocities  at  the  largest  ranges  where  void  crushup  is  incomplete. 

2.  Shock  damaged  shear  modulus:  A  decrease  in  shear 
modulus  of  a  factor  of  almost  three  during  loading  resulted  in  larger 
peak  particle  velocities  (in  worse  agreement  with  the  measurements), 
slightly  wider  positive  velocity  pulses,  and  approximately  50  percent 
wider  negative  velocity  pulses. 

3.  Shock  damaged  failure  surface:  The  failure  surface  used 
for  previous  work  was  based  on  measurements  of  stress  differences  at 
confining  pressures  of  4  Kb.  A  calculation  was  made  in  which  the 
allowable  stress  difference  after  shock  loading  was  decreased  40%  . 
This  was  based  on  laboratory  material  properties  measurements  made  at 
TerraTek.  Strain  rate  dependence  was  assumed  to  Increase  the  strength 
of  damaged  grout  In  a  similar  manner  to  undisturbed  grout.  The  shock 
damaged  failure  surface  resulted  in  positive  particle  velocity  pulses 
wider  than  the  measurements,  peak  velocities  larger  than  the  data,  an 
18%  reduction  In  calculated  peak  residual  hoop  stress  and  a  slightly 
larger  cavity.  Hot  much  Improvement  in  negative  velocity  pulse  width 
was  achieved.  On  the  whole  however,  the  concept  of  a  damaged  failure 
surface  appeared  to  be  an  Improvement  in  the  constitutive  model. 
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4.  P-V  EOS  and  elastic  moduli  of  grout:  These  variations 
gave  second  order  changes  In  results  at  most. 

5.  Strength,  strain  rate  dependence,  and  EOS  changes  In  the 
models  for  the  luclte  shell  surrounding  explosive:  These  variations 
gave  only  second  order  changes  In  results. 

6.  HE -detonated  or  burned:  At  ranges  greater  than  1  cm, 
almost  no  changes  In  results  from  the  full  HE  detonation  calculation 
occurred  when  the  HE  was  assumed  to  burn  at  constant  volume  and  then  to 
follow  the  JWL  equation  of  state  Isentrope  down  In  pressure  from  the 
Initial  volume. 

7.  Modified  JWL  constants  for  PETN  explosive:  Changes  in  the 
equation  of  state  of  the  PETN  explosive  (density  3  1.0  gm/cc)  that  did 
not  alter  the  measured  Chapman-Jouget  pressure,  detonation  velocity  and 
explosive  energy  were  Investigated.  Some  reduction  in  peak  particle 
velocities  was  achieved  at  the  velocity  gauge  locations  Inside  2.0  cm. 


Based  on  the  above  constitutive  modeling  variations,  particu¬ 
larly  those  which  improved  the  calculatlonal  results  significantly, 
l.e.,  changes  (2)  and  (3),  we  hypothesized  that  a  weakened  central  core 
of  grout  was  needed  to  give  the  desired  wide  negative  velocity  pulses. 
All  of  the  following  modeling  variations  were  directed  toward  that 
objective. 

8.  No  strain  rate  strength  dependence  for  damaged  grout: 

This  further  reduction  In  the  strength  of  the  shock  damaged  grout  (from 
(3))  resulted  In  higher  peak  velocities  and  wider  positive  velocity 
pulses,  both  of  which  were  unsatisfactory,  but  widened  the  negative 
pulses  significantly. 


9.  Other  variations  which  either  enlarged  the  damaged  central 
core  or  Increased  the  amount  of  damage  beyond  that  measured  In  the 
laboratory  widened  the  negative  velocity  pulses  significantly  but  also 
widened  the  positive  pulses  as  well. 

10.  Plastic  work  effects;  The  above  results  Indicated  that 
the  grout  remained  strong  during  shock  loading  and  only  decreased  In 
strength  later.  If  true,  this  would  result  in  a  sufficiently  narrow 
positive  velocity  pulse  and  a  wide  negative  pulse.  As  a  first  try,  the 
strain  rate  dependence  on  strength  was  made  a  function  of  the  magnitude 
of  the  plastic  work  in  a  grout  element,  more  plastic  work  resulting  In  a 
gradually  lower  dependence  of  strength  upon  strain  rate.  At  some 
prescribed  value  of  plastic  work,  no  strain  rate  dependence  was  allowed 
for  the  damaged  grout.  This  gave  a  50  percent  increase  In  negative  pulse 
width  compared  to  the  results  with  standard  strain  rate  dependence  and 
shock  damage. 

11.  Time  dependent  damage:  In  this  model,  the  undamaged 
(Intact  grout)  failure  surface  was  assumed  to  be  strain  rate  dependent 
while  the  damaged  failure  surface  was  given  no  strain  rate  dependence. 

A  time  dependent  (Maxwell  solid)  relaxation  of  devlatoric  stresses  from 
the  undamaged  to  the  damaged  failure  surface  was  assumed  to  occur  In  a 
relaxation  time  dependent  upon  the  amount  of  damage.  This  model  has 
given  the  best  results  to  date,  l.e.,  correct  peak  velocities  and 
positive  velocity  pulses,  and  negative  velocity  pulses  only  slightly 
narrower  than  the  measurements.  At  years  end,  this  constitutive  model 
was  being  pursued  further  with  excellent  prospects  for  a  successful 
simulation  of  the  velocity  data. 

5.2  SANDIA  HE  TESTS  IN  G-TUNNEL  TUFF 

Sandla  National  Laboraborles  have  been  conducting  high  explo¬ 
sive  tests  In  tunnel  bed  5  of  G-tunnel,  Area  12,  at  NTS.  Radial 
stresses  at  several  ranges  from  the  explosive  and  cavity  pressure  have 
been  measured  out  to  times  of  many  hours  as  an  experimental  test  of  the 
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residual  stress  concept.  The  largest  of  these  tests.  Puff  Too,  described 
by  Smith,  Reference  (7),  used  1000  1b  of  TNT  explosive.  Stress  vs  time 
data  for  comparison  with  dynamic  calculations  was  obtained  from  Puff  Too 
as  well  as  from  two  smaller  64  1b  tests,  RS11  and  RS15.  We  have  conducted 
finite  difference  simulations  of  these  tests  to  Improve  our  understanding 
of  constitutive  models  for  tuff. 

The  Sandia  stress  measurements  showed  a  considerable  scatter  from 
test  to  test  as  expected  for  an  earth  material  as  variable  as  the  satura¬ 
ted  tuff.  However,  several  features  of  the  experimental  results  were 
quite  Interesting.  The  1000  1b  event  gave  an  average  cavity  radius  of 
108  cm,  which  could  easily  be  calculated  with  our  constitutive  models. 
However,  the  64  1b  tests  gave  a  significantly  smaller  cavity  (approxi¬ 
mately  82  cm  when  yield  scaled  to  1000  1b).  Calculations  showed  that 
more  than  65%  of  the  difference  In  scaled  cavity  radius  could  be 
explained  by  using  the  same  failure  strength  for  both  If  a  strain  rate 
dependent  strength  model  was  Invoked.  Strain  rate  dependent  strength 
Increases  had  already  been  used  to  explain  the  small  cavity  radii 
measured  In  SRI  laboratory  tests  In  grout  and  In  tuff.  A  smaller  scale 
test  would  have  higher  strain  rates  than  a  larger  test  and  therefore  a 
higher  strength  and  smaller  cavity. 

Another  feature  of  Interest  was  the  width  of  a  dip  In  the  stress 
vs  time  measurements  seen  both  In  RS11  and  RS15.  This  dip  has  been 
correlated  through  numerical  simulations  with  the  width  of  the  negative 
velocity  pulse,  l.e.,  with  the  duration  of  the  cavity  rebound.  A 
calculation  was  made  Introducing  a  damage  dependent  decrease  in  shear 
modulus  to  widen  the  negative  velocity  pulse  which  also  improved  the 
calculated  width  of  this  dip  In  the  stress  record.  When  the  velocity 
records  from  ther  SRI  grout  spheres  tests  are  successfully  simulated, 
calculations  of  the  Sandia  tuff  tests  made  with  these  constitutive 
models  can  be  expected  to  widen  the  stress  dip  as  well. 


5.3  RESIDUAL  STRESS  RELAXATION 


A  central  concern  of  the  DNA  containment  program  over  the  last 
year  or  so  has  been  to  build  on  the  earlier  work  which  demonstrated 
the  existence  of  a  residual  stress  field  and  determined  how  that  field 
may  change  In  time.  The  SRI  work  has  clearly  demonstrated  that  if 
hydrofracture  of  an  exploded  grout  sphere  Is  delayed  for  a  period  of 
many  minutes,  the  measured  hydrofracture  pressure  becomes  comparable 
to  that  for  an  unexploded  sphere.  In  other  words,  a  residual  stress 
has  relaxed  away  by  some  process  or  other.  A  similar  conclusion  can 
be  drawn  from  hydrofrac  experiments  at  NTS  where  no  evidence  of  a 
residual  stress  field  around  a  former  nuclear  test  can  be  found  some 
months  after  the  experiment. 

Two  mechanisms  for  stress  relaxation  have  been  suggested;  one 
involves  material  creep,  the  other  pore  fluid  migration.  We  made  a 
number  of  calculations  several  years  ago  with  a  very  simple  effective 
stress  law  which  assumed  that  the  pore  pressure  was  equal  to  the  mean 
stress  for  fully  saturated  materials.  Under  these  circumstances,  the 
residual  stress  was  calculated  to  relax  away  on  a  time  scale  governed 
by  the  permeability  of  the  material  and  the  viscosity  of  water.  How¬ 
ever,  more  elaborate  effective  stress  models  have  been  developed  for 
elastic  materials.  We  have  attemped  to  incorporate  these  represen¬ 
tations  In  order  to  do  a  "better"  job  of  representing  the  physical 
situation. 

In  order  to  make  progress,  It  Is  necessary  to  determine  a 
number  of  material  moduli  which  deal  with  drained  and  undrained 
materials.  These  moduli  are  not  determined  in  the  conventional  suite 
of  measurements  made  at  TerraTek,  so  an  experimental  program  was  under¬ 
taken  to  determine  the  needed  values.  Unfortunately  the  time  scale 
for  the  attainment  of  dralned-materlal  equilibrium  was  found  to  be 
very  long,  and  there  Is  serious  question  about  the  reliability  of  the 
experimental  determinations.  Nevertheless,  a  series  of  stress  relaxa¬ 
tion  calculations  was  made  with  the  best  estimates  of  material  moduli 
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for  2C4  grout.  As  a  result  the  calculated  residual  stress  relaxation 
would  reduce  the  peak  compressive  stress  around  an  explosive  cavity  by 
only  about  15%  .  This  value  Is  disappointingly  small  compared  to  the 
100%  relaxation  determined  by  the  simpler  model  and  apparently  required 
by  experiment. 

At  the  present  time  further  Investigation  Is  underway  to  determine 
the  best  effective  stress  law  treatment.  The  relatively  simple  model 
previously  used  can  be  questioned  for  elastic  materials.  The  full 
elastic  treatment  can  also  be  questioned  because  there  Is  reason  to 
believe  that  material  damage  which  would  be  particularly  severe  near 
the  explosive  cavity  would  significantly  modify  material  properties  and 
bring  Into  question  the  applicability  of  an  elastic  treatment  for  that 
material.  More  definitive  results  are  anticipated  within  the  next 
reporting  period. 

5.4  MISCELLANEOUS  ACTIVITIES 

We  have  maintained  a  very  active  Interest  In  many  aspects  of  the 
grout  spheres  program  through  the  years.  As  a  result  a  number  of 
analyses  of  the  available  data  have  been  made,  and  suggestions  for 
additional  work  have  been  made.  These  will  be  summarized  very  briefly 
below. 


There  was  some  question  concerning  the  reliability  of 
the  first  set  of  decoupled  tests  made  last  year.  Data  from 
these  experiments  was  carefully  analyzed,  and  It  appeared 
very  likely  that  the  precast  cavities  were  filled  with 
water  before  the  explosive  charge  was  fired.  The  experi¬ 
mental  technique  has  been  modified  and  the  experiments 
redone  successfully. 


An  overall  evaluation  of  the  contribution  of  this  work  has 
been  prepared  for  the  COR  in  which  the  demonstration  of  the 
existence  of  the  residual  stress  field  and  Its  relaxation 
in  rock- type  materials  has  been  emphasized. 
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•  We  have  contributed  to  periodic  reviews  of  the  grout  sphere 
program. 

•  An  experiment  was  suggested  In  which  the  explosive  charge 
would  be  fired  In  a  small,  closed  metal  vessel.  This  was 
proposed  as  a  technique  to  assure  that  all  aspects  of  the 
experimental  apparatus  performed  as  expected.  This 
experiment  has  been  conducted  successfully. 

•  Since  two  different  mechanisms  have  been  suggested  to 
explain  the  observed  residual  stress  relaxation.  It  seemed 
appropriate  to  suggest  experiments  which  could  choose 
between  the  mechanisms.  The  most  straightforward  possi¬ 
bility  would  Involve  experiments  on  dry  material  where 
pore  fluid  diffusion  Is  not  Involved.  To  this  end 
suggestions  were  made  which  might  be  suitable  for  the 
preparation  of  a  dry  alluvium- like  material.  In  addition, 
other  experiments  were  devised  which  could  have  shown  the 
Influence  of  creep  and  pore  fluid  diffusion,  but  the  one 
dimensional  beam  geometry  suggested  was  far  from  that  of 
Immediate  concern  In  the  containment  program. 
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SECTION  6 
MISCELLANEOUS 


A  number  of  tasks  were  undertaken  during  the  year  which  do  not 
fall  well  within  the  topics  discussed  above.  They  will  be  grouped 
together  here. 

6.1  KRAK  MODELING  ASSISTANCE 

The  KRAK  code  has  been  used  at  LANL  to  Investigate  fracture 
propagation  from  a  borehole.  Because  the  Sneddon  crack  shape  formalism 
Is  used  by  that  code,  some  approximation  must  be  made  to  the  pressure 
distribution  from  the  origin  to  the  cylindrical  cavity  wall  In  order 
to  perform  the  Sneddon  numerical  Integration.  It  Is  desirable  to 
optimize  the  pressure  distribution  assumption  so  that  the  crack  shape 
calculated  by  the  Sneddon  Integral  agrees  as  closely  as  possible  with 
the  "correct"  shape  as  determined  by  finite  element  calculations  for 
the  real  geometry  being  considered  and  making  use  of  the  pressure 
distribution  thought  to  exist  on  crack  surfaces.  The  following  discus¬ 
sion  describes  the  work  done.  First  Is  a  treatment  of  the  analytical 
and  numerical  solution  to  the  Sneddon  problem.  This  Is  followed  by  a 
discussion  of  the  finite  element  calculations  of  the  expected  crack 
shape  In  the  vicinity  of  a  pressurized  borehole  and  finally,  the 
recommended  numerical  algorithm  to  be  used  In  KRAK  to  reproduce  the 
finite  element  solution. 

6.1.1.  Analytical  Crack  Nldth  Calculation 

As  part  of  the  effort  to  provide  an  Improved  algorithm  for  crack 
width  calculations  In  KRAK,  several  analytical  crack  width  calculations 
were  done  using  the  LANL  pressure  distribution  Inside  the  crack  and 
the  Sneddon  formalism.  In  an  Infinite  elastic  medium  the  crack  width 
can  be  represented  by  Integrals  of  the  pressure  distribution  along  the 
crack  wall. 


Rfi 


Four  different  cases  were  considered: 


CASE  I:  Linear  crack  of  length  L.  This  Is  expected  to  be  a 

good  approximation  when  the  length  of  crack  Is  short  compared  to  the 
cavity  size.  Assuming  plane  strain  conditions,  the  half  width  of  the 
crack  at  the  position  x  along  the  crack  Is  given  by 


In  the  above  equation  G  and  v  represent  the  shear  modulus  and  Poisson's 
ratio  of  the  material  and  p(x)  represents  the  pressure  distribution 
along  the  crack. 


CASE  II:  Penny  shaped  crack  with  the  radius  equal  to  the 

crack  length  L.  The  half  crack  width  at  the  distance  r  from  the  axis 
Is  given  by 


W(r) 
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where  R  ■  L  and  fp  ■  ^  . 


CASE  III:  Large  penny  shaped  crack  with  the  cavity  region  at 

the  origin  without  cavity  pressure.  The  half  crack  width  Is  given  by 
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In  the  above  equations,  R  represents  the  cavity  radius. 

CASE  IV;  Same  as  CASE  III  with  the  cavity  pressure  applied 
In  the  cavity  region.  The  Integral  expression  for  CASE  II  can  be 
directly  used  with  R  *  R£  +  L. 

For  the  numerical  Integration  the  pressure  was  assumed  to  vary 

linearly  within  each  cell,  l.e..  In  the  1-th  cell  with  p.  and 

/rl  M 

fj  *  I  jp  or  given  at  the  left  and  p^+1  and  f1+l  given  at  the 
right  boundary  of  cell, 

p(f)  *  +  b|  f  for  ff*f*fj+1 

Pl+1  *  Pf 

where  '  Pj  -  fm  -''fj  fi 
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Using  this  approximation  the  evaluation  of  the  first  Integrals  at  the 

J-th  point  becomes  the  sum  of  the  following  Integrals 


r1 

■i  J. 


1+1  df  f,+1  f  df 

,  Vf/  -  7  ’  1  f,  ^ 


and  y 


f2df 

1  VfjZ  -  f‘ 


from  1  ■  I  to  J  -  1 


60 


\ 

- 4 


I 


ti-wuM 

S 


J 

Jra* 


A  :A 


US* 


•'a' 


*1 


••Aja 


where  a^'s,  6 ^ ' s*  and  * s  are  constant  within  each  cell.  These 
elementary  Integrals  can  be  easily  evaluated.  After  the  first  of  the 
double  Integrals  Is  evaluated  at  each  discrete  point*  these  Integrals 
are  again  assumed  to  vary  linearly  within  each  cell.  Following  the 
same  procedure  the  second  Integral  at  the  K-th  point  can  be 
approximated  by  the  sum  of  the  following  Integrals, 


fra.  1  *  K  to  IMX  -  1 

where  a^'s,  and  r^'s  are  constant  within  each  cell. 

The  results  of  crack  width  calculations  using  the  LANL  pressure 
distribution  for  the  four  cases  are  shown  In  Figure  6.1.  It  is  clear 
that  the  Initial  LANL  result  corresponded  to  Case  II  which  Ignored  the 
presence  of  the  borehole  and  assumed  the  crack  simply  translated  to 
the  axis.  It  Is  also  clear  that  other  assunptlons  about  how  to  solve 
the  Sneddon  problem  In  this  real-world  case  lead  to  very  different 
estimates  of  crack  width  for  a  given  pressure  distribution.  However, 
the  crack  width  Influences  flow  into  the  crack,  the  crack  length,  and 
the  pressure  distribution,  so  It  Is  Important  to  use  the  best  approxi¬ 
mation  available  for  the  crack  width  formalism. 

6.1.2.  Simple  Algorithm  for  Crack  Width  Calculation 

The  next  step  was  to  determine  the  "correct"  crack  width  for 
the  actual  geometry  under  Investigation  by  making  finite  element  calcu¬ 
lations.  This  provided  a  standard  of  comparison  for  the  algorithms  to 
be  developed.  Figure  6.2  shows  the  crack  width  calculated  for  a  5.9 
cm  long  crack  extending  out  from  a  cylindrical  cavity  of  7.5  cm  radius 
and  60  cm  length.  The  LANl-supplled  pressure  distribution  was  used  In 
the  crack. 
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Figure  6.1.  Crack  width  comparison  for  LASL  Pressure  Distribution. 
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6.2.  Comparison  of  Calculated  Crack  Widths 
for  Different  Pressure  Distributions  In  Cavity 


Effective  Pressure 


For  comparison,  the  several  different  pressure  distributions 
shown  In  Figure  6.2  were  also  used  In  Sneddon  calculations  to  determine 
the  various  shapes  shown  In  the  figure.  The  simple  expedient  of 
extrapolating  the  pressure  profile  0.6  cm  Into  the  borehole  (c),  or 
extending  the  cavity  pressure  1  cm  Into  the  borehole  (d)  give  good 
agreement  with  the  correct  crack  shape  within  the  uncertainty  of  the 
finite  element  method. 

The  question  next  arises  as  to  how  the  approximate  algorithms 
work  for  longer  cracks.  This  was  Investigated  by  simply  stretching 
out  the  LANL  pressure  distribution  to  crack  lengths  of  29.5  and  59 
cm.  Figure  6.3  shows  a  comparison  of  several  analytic  crack  width 
results  with  the  "correct"  finite  element  solution.  In  this  and  the 
following  figures  the  "pressure  fix"  referred  to  In  the  legend  means 
that  the  cavity  pressure  was  assumed  to  act  on  1  cm  of  the  borehole  at 
the  cavity  wall  as  In  profile  d  of  Figure  6.2.  Figure  6.4  shows 
similar  results  for  a  59  cm  long  crack. 

Figure  6.5  demonstrates  the  relative  Importance  of  cavity 
Influence  for  the  different  crack  lengths.  The  ratio  of  the  crack 
width  at  the  opening  for  the  case  of  cavity  pressure  applied  to  that 
for  the  case  without  cavity  pressure  applied  and  the  ratio  of  the 
crack  opening  of  real  case  to  that  of  penny  shaped  crack  without 
cavity  pressure  are  shown  for  the  three  different  crack  lengths.  The 
ratio  of  the  real  crack  width  to  the  width  calculated  with  the  cavity 
pressure  In  1  cm  of  the  cavity  Is  also  shown  In  this  figure. 

These  results  demonstrate  that  the  simple  scheme  of  loading  the 
cavity  pressure  In  1  cm  of  the  crack  adjacent  to  the  opening  of  a 
penny  shaped  crack  leads  to  analytic  results  which  approximate  the 
correct  shape  quite  well  even  for  cracks  5  to  10  times  longer  than 
those  actually  calculated  by  LANL.  This  suggests  a  simple  calcula- 
tlonal  procedure  to  be  used  In  KRAK  to  calculate  crack  openings  from 
this  particular  cylindrical  borehole  geometry. 


64 


figure  6.3  Comparison  of  Calculated  Crack  Widths 

for  Different  Geometry  Cases  (29.5  cm.  crack ) 


Radial  direction  in  centimeters 


figure  6.4.  Comparison  of  Calculated  Crack  Widths 
for  Different  Geometry  Cases  (59.  cm.  era 
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Radial  direction  in  centimeters 


figure  6.5.  Cavity  Influence  as  a  Function  of  Crack  Opening 


Radius 


This  calculated  algorithm  was  developed  for  this  particular 
geometry.  There  Is  no  reason  to  expect  It  to  apply  to  any  other 
problem  such  as  a  crack  propagating  from  a  spherical  cavity*  The 
whole  process  of  comparing  finite  element  results  and  analytic 
approximations  Including  various  treatments  of  the  cavity  region  must 
be  repeated  in  order  to  find  an  appropriate  algorithm  for  a  different 
problem. 

6.2  CONTAINMENT  SYMPOSIUM 

S-CUBED  personnel  presented  three  papers  at  the  Containment 
Symposium  held  at  the  Naval  Post-Graduate  school,  Monterey, 

California,  August  26-28,  1981.  These  papers  are: 

"Medium  Property  Effects  on  Containment"  -  Norton  Rimer 
"The  Nature  of  LOS  Pipe  Flow"  -  J.  R.  Barthel 
"Evaluation  of  Red  Hot"  -  Russell  E.  Duff 

These  three  documents  will  appear  in  the  Proceedings  of  the  Symposium. 

6.3  CEP  PARTICIPATION 

A  member  of  the  S-CUBED  Staff  has  continued  to  serve  as  a 
consultant  member  to  the  Department  of  Energy  Containment  Evaluation 
Panel.  This  organization  reviews  and  comments  on  the  plans  for  all 
underground  nuclear  tests  to  be  conducted  by  the  United  States.  This 
Involvement  has  represented  a  relatively  minor  but  continuing  expendi¬ 
ture  of  contract  resources. 
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